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Monitoring resistance to Bacillus thuringiensis subsp. israelensis in 
the field by performing bioassays with each Cry toxin separately 

Guillaume Tetreau/% Renaud Stalinski, Jean-Philippe David, Laurence Despres 

Laboratoire d'Ecologie Alpine, Unite Mixte de Recherche 5553, Universlte de Grenoble 1, Grenoble, France 

Bacillus thuringiensis subsp. israelensis (^Bti) is increasingly used worldwide for mosquito control and is the only 
larvicide used in the French Rhone-Alpes region since decades. The artificial selection of mosquitoes with field- 
persistent Bti collected in breeding sites from this region led to a moderate level of resistance to Bti, but to relatively 
high levels of resistance to individual Bti Cry toxins. Based on this observation, we developed a bioassay procedure 
using each Bti Cry toxin separately to detect cryptic "Qii-resistance evolving in field mosquito populations. Although 
no resistance to Bti was detected in none of the three mosquito species tested fAedes rusticus, Aedes sticticus and 
Aedes vexansj, an increased tolerance to Cry4Aa (3. 5 -fold) and CryllAa toxins (8-fold) was found in one Ae. sticti- 
cus population compared to other populations of the same species, suggesting that resistance to Bti may be arising 
in this population. This study confirms previous works showing a lack o/Bti resistance in field mosquito popula- 
tions treated for decades with this bioinsecticide. It also provides a first panorama of their susceptibility status to 
individual Bti Cry toxins. In combination with bioassays with Bti, bioassays with separate Cry toxins allow a more 
sensitive monitoring ofBti-resistance in the field. 
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Numerous insect species represent major econom- 
ic and sanitary issues for human societies as they can 
strongly impact crop cultures or transmit severe human 
diseases such as malaria, leishmaniasis or trypano- 
somiasis (Medlock et al. 2012, Savary et al. 2012, van 
Griensven & Diro 2012). The control of pest populations 
is therefore of high importance and relies mainly on the 
use of chemical and biological insecticides (Meissle et al. 
2011, van den Berg 2011). The main threat for an efficient 
pest control is the ability of insects to quickly develop 
resistance to nearly all the insecticides used, leading to a 
strong decrease of insecticides efficiency (Tabashnik et 
al. 2008, Rivero et al. 2010). To bypass these resistances 
and adapt the insecticide treatment strategies, detecting 
incipient resistance in treated populations as early as 
possible is required (Carriere et al. 2010). 

The most obvious approach to detect resistance in 
the environment is to perform bioassays on field-sam- 
pled individuals, which consist in exposing insects to 
increasing amount of insecticide to calculate their rela- 
tive tolerance to this insecticide (Sivasupramaniam et 
al. 2007). Such basic approach is often preliminary to 
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further studies exploring in-depth the molecular bases 
of the resistance (Tabashnik et al. 2006, Marcombe et 
al. 2009). Bioassays have been performed on numerous 
field populations of insects with different chemical in- 
secticides (Ocampo et al. 2011, Shin et al. 2011, Silva 
et al. 2011, Springate & Colvin 2012), but also with bio- 
logical insecticides such as Bacillus sphaericus (Su & 
Mulla 2004, Silva et al. 2008, Akiner et al. 2009) and 
Bacillus thuringiensis subspecies (Tabashnik et al. 1990, 
Gao et al. 2010, Tabashnik & Carriere 2010, Alcantara et 
al. 2011, Gassmann et al. 2011). B. thuringiensis subsp. 
israelensis (Bti) is increasingly used worldwide for mos- 
quito control and investigating potential Bti-resistance 
in the environment is essential to ensure the long term 
use of this insecticide. Resistance to Bti has been ex- 
tensively searched in natural populations of mosquitoes 
(Liu et al. 2004, Vasquez et al. 2009, Loke et al. 2010, 
Kamgang et al. 2011) and generally concluded in a lack 
of Bti resistance in field populations even after decades 
of Bti treatments, although a few studies described a de- 
creased susceptibility in some populations (Zhang et al. 
2004, Paul et al. 2005, Boyer et al. 2007, 2012). Never- 
theless, these studies may suffer from technical biases, 
notably on the choice of the reference strain. 

Usually, when bioassays are performed, a laboratory 
strain from the same species, susceptible to the insec- 
ticide tested, is used as a reference to calculate the re- 
sistance ratio of the strains/populations tested. Although 
this approach is widely used, it is adapted only if the ref- 
erence strain has been sampled recently in populations 
geographically close to those tested. If this is not the 
case, results can reflect a basal difference of tolerance 
between the field populations and the reference strain, 
independent from fi?/-treatments, rather than a selected 
resistance evolving in the populations treated with Bti 
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(Zhang et al. 2004, Paul et al. 2005). Moreover, this ap- 
proach is only possible for a restricted number of species, 
as most mosquito species cannot be maintained in labo- 
ratory conditions. An alternative solution is to compare 
treated and untreated populations; but even in this case, 
the differences observed can be due to natural tolerance 
differences between populations and not to a selected 
resistance (Boyer et al. 2007, 2012). Therefore when 
mosquitoes cannot be maintained in laboratory, the only 
way of detecting incipient resistance is to apply a highly 
sensitive bioassay approach on field-populations and test 
for the evolution of their tolerance for several years (Jan- 
maat & Myers 2003, Tabashnik & Carriere 2010). Strong 
evidence for an evolved resistance in the environment is 
given by an increasing tolerance to the insecticide in the 
same population over generations. Unfortunately, one of 
the main reasons why no firm resistance to Bti was de- 
tected in the field yet may be due to the low resolution 
of bioassays performed with Bti, a mixture of four main 
toxins (Cry4Aa, Cry4Ba, CryllAa and CytlAa). 

Recent works on Bti-resistance showed that while 
levels of resistance to Bti can be low in a laboratory-se- 
lected strain (3.5-fold), resistance to each of the Cry tox- 
ins tested separately in bioassays can be much higher (up 
to 60-fold) (Paris et al. 2011, Tetreau et al. 2012b). Here 
we present an improved bioassay method consisting in 
performing distinct bioassays with each of the three Bti 
Cry toxins on field populations of mosquitoes. We ap- 
plied this approach to various mosquito populations col- 
lected in the French Rhone-Alpes region, where Bti is the 
only larvicide used for mosquito control since 1990. Fur- 
thermore, the leaf litters containing field-persistent Bti 
able to select for resistance in mosquitoes in the labora- 
tory were sampled in this region, raising concerns about 
possible 5?/-resistance evolving in mosquitoes in these 
particular breeding sites (Tilquin et al. 2008, Paris et 
al. 2010, 2011). The variations in susceptibilities of field 
mosquito populations to Bti and Cry toxins are described 
and discussed. Perspectives of further improvement and 
wide application of this tool are also discussed. 

MATERIALS AND METHODS 

Mosquito sampling and species identification - Mos- 
quito larvae were sampled in 17 different breeding sites 
of the French Rhone-Alpes region in close collaboration 
with the Entente Interdepartementale pour la Demousti- 
cation (EID) (Fig. 1). Among them, 12 sites were regu- 
larly treated with Bti for decades and five had never been 
treated with Bti. Each sampling site is identified by the 
letter "P" followed by a number. Each site can be associ- 
ated to different "P" numbers corresponding to different 
sampling dates. As mosquito breeding sites often con- 
tained more than one mosquito species, different species 
sampled at the same time in the same site were associ- 
ated to the same "P" number. 

As numerous different species can share the same 
breeding site, larvae were individually identified before 
bioassays using the identification key provided in Becker 
et al. (2010). Bioassays with different species were per- 
formed independently. 



Production of individual Bti Cry toxins - A crystal 
negative strain of Bti (4Q2-81) transformed with the 
plasmids pHT606, pHT618 or pWF53 was used for the 
production of Cry4Aa, Cry4Ba or CryllAa toxins re- 
spectively. These strains were obtained from the Pasteur 
Institute (Paris, France) or from Prof B Federici (Uni- 
versity of California, Riverside, CA, USA). Each toxin 
was produced, purified and quantified as previously de- 
scribed (Tetreau et al. 2012b). Toxin suspensions were 
stored at -20°C until use. 

Bioassays procedure and data analyses - To avoid 
any influence of the larval stage on bioassays, larvae 
were maintained in the water from the breeding site 
where they were sampled until the third instar. Consid- 
ering that the species studied cannot reproduce in labo- 
ratory conditions, we could not perform bioassays on Fl 
larvae, but directly used field-collected larvae. In order 
to limit the impact of environmental parameters in the 
sampled sites {e.g. pH of the water, pollutants etc.) on 
bioassay results, field collected larvae were kept in fresh 
water for 24 h before performing the bioassays. 

When more than 2,000 larvae of the same species 
were sampled, standard bioassays were conducted by 
exposing larvae to five-seven concentrations of Cry4Aa, 
Cry4Ba, CryllAa and commercial Bti (Vectobac WG, 
3500 ITU/mg) for 48 h in order to obtain 5-95% mortal- 
ity. When 1,000-2,000 larvae were sampled, high sensi- 
tive standard bioassays could be performed, but not for 
all the Cry toxins; bioassays with Bti and Cry4Aa were 
performed in priority since Cry4Aa is the main toxin 
persisting in the field after Bti treatment (Tetreau et al. 
2012a) and with the highest resistance ratio in the Bti- 




Fig. 1: localisation of the 32 mosquito breeding sites visited in the 
French Rhone-Alpes region. The 22 treated (triangles) and the 10 
untreated (circles) sites visited are shown. Sampled sites used for 
bioassays (12 for treated and 5 for untreated sites) are indicated by a 
filled item and by the name of the population in italic (untreated site) 
or not (treated site). 
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resistant strain described in Tetreau et al. (2012b). Lethal 
concentrations for 50% individuals (LC^j,) were calcu- 
lated for each population and each toxin using a probit 
statistical model with the module "dose" of XLSTAT 
v.2009.4.06 (Addinsoft). 

When only a limited number of larvae of the same 
species were sampled, standard bioassays could not be 
performed, as a reduced number of replicates and/or con- 
centrations performed would lead to a high variability in 
larval mortality and to a low reliability of the LC^^ val- 
ues. In this case, a "diagnostic dose" approach was used, 
consisting in using one concentration of Cry toxin or Bti 
chosen in order to obtain 60-80% mortality. These doses 
were designed for each species and each toxin from the 
bioassays data from the first population sampled. For 
each toxin and Bti, statistical differences in larval mor- 
tality between two populations from the same species 
were measured by a Wilcoxon test performed with R 
2.8.1 software (R Development Core Team 2007). 

For each toxin tested, four replicates per concentra- 
tion was performed on 20 third-instar larvae in 50 mL 
of tap water containing the insecticide according to the 
standard bioassay procedure described by the World 
Health Organization (WHO 2005). For each population 
tested, five negative controls, consisting in 20 larvae ex- 
posed to tap water only, were performed. All bioassays 
were performed at 17°C. For each population, the mean 
larval size was calculated from more than 15 larvae. 

Significance of larval size and mortality differences 
between the three species (Aedes rusticus, Aedes sticti- 
ciis and Aedes vexans) were estimated by performing a 
Wilcoxon test using R 2.8.1 software (R Development 
Core Team 2007). 

Considering that the species studied cannot be main- 
tained in laboratory, we were unable to compare the tol- 
erance levels of the sampled populations to a susceptible 
laboratory strain. Therefore, for each species we com- 



pared populations of mosquitoes with each other. Then, 
we calculated "tolerance ratios" which represent an in- 
creased or decreased susceptibility of a population rela- 
tive to the others. 

RESULTS AND DISCUSSION 

Limits of performing bioassays on field populations 
- Among the 71 natural mosquito populations sampled, 
data from 59 populations (83%)) were not analysable due 
to insufficient larvae to perform bioassays (69%)), to mor- 
tality in controls (11%)) or to high rate of pupation dur- 
ing bioassays (3%o) (Table I). No bioassay data could be 
obtained for Aedes cantons, Aedes cataphylla, Aedes ci- 
nereus and Culex pipiens field populations in the present 
study, while sufficient data were obtained forAe. rusticus, 
Ae. sticticus andAe. vexans populations (Table I). Among 
the three Ae. rusticus populations tested in bioassays, two 
were sampled in treated sites (P3 and P5) and one in an 
untreated site (P14) while ^e. sticticus (P16, P20 and P52) 
and^e. vexans (P18, P20, P30, P42 and P44) populations 
were all sampled in treated sites. Only one population of 
Culiseta annulata (P51) was sampled in sufficient amount 
to perform a diagnostic dose approach and showed com- 
parable mortality rates as Ae. vexans populations at the 
same doses (data not shown). 

The fact that only less than 17% of the populations 
sampled and only four out of the eight species sampled 
gave analysable results highlights the difficulty to find 
sites producing enough calibrated larvae of the same 
species to perform bioassays. This limit is mainly due 
to typical features of the Rhone -Alpes region, consti- 
tuted of numerous sites of moderate size largely scat- 
tered in woodlands. 

Interspecific comparisons - Ae. sticticus and Ae. 
vexans larvae exhibited the same patterns of tolerance 
(Wilcoxon test, p > 0.05) except for Cry4Aa for which 



TABLE I 

Summary of mosquito breeding sites sampled where at least one of the eight studied species were found 



Bioassays Not analysable populations 



Species 


Sampled 




Diagnostic 
dose 


Low 
sample size 


High mortality 
in controls 


Pupation 


Aedes cantans 


11 


0 


0 


9 


2 


0 


Aedes cataphylla 


5 


0 


0 


5 


0 


0 


Aedes cinereus 


3 


0 


0 


1 


1 


1 


Aedes rusticus 


15 


3 


0 


12 


0 


0 


Aedes sticticus 


14 


3 


0 


8 


3 


0 


Aedes vexans 


12 


3 


2 


6 


0 


1 


Culex pipiens 


6 


0 


0 


4 


2 


0 


Culiseta annulata 


5 


0 


1 


4 


0 


0 


Total 


71 


9 


3 


49 


8 


2 



for each species, the number of populations for which bioassays were performed is indicated. The main reason why bioassays 
were not performed for the other populations is also indicated. LC^u! lethal concentrations for 50% individuals. 
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TABLE II 

Mean lethal concentrations for 50% individuals (LCj^) calculated for all the populations of each species and tolerance ratios 
obtained for each species with individual Cry toxins and commercial Bacillus thuringiensis subsp. israelensis {Bti) 

Size Cry4Aa Cry4Ba CryllAa Bti" 

Species (mm) (ng/mL) (ng/mL) (ng/mL) (|ig/mL) 

Aedesnisticus 10.6 ±0.8 262.7 ±88.3 60.9 ±15.2 106.9 ±22.2 2.03 ±1.05 

Aedes sticticus 6.3 ± 0.9 14.8 ±6.8 8.9 ± 0.0 15.7 ±11.0 0.11 ±0.028 

Aedesvexans 6.2 ± 0.6 67.6 ± 6.4 13.5 ±0.4 27.1 ± 11.7 0.11 ±0.009 

Tolerance ratios 

Ratio spl/sp2 Size Cry4Aa Cry4Ba CryllAa Bti" 

Ae. rusticus/Ae. sticticus 1.7x 17.8x 6.8x 6.8x 18.5x 

Ac. rusticus/Ae.vexans 1.7x 3.9x 4.5x 3.9x 18. 5x 

Ae. vexans/Ae. sticticus l.Ox 4.6x 1.5x 1.7x l.Ox 

a: Vectobac WG. Larval size ratios between species are also indicated. 



Ae. vexans larvae were 4.6-foId more tolerant than Ae. 
sticticus larvae (Wilcoxon test, p < 0.05) (Table II). Ae. 
rusticus was significantly 3. 9-fold-18. 5-fold more toler- 
ant than^e. sticticus andAe. vexans species when tested 
with Cry toxins and commercial Bti (Wilcoxon test, p < 
0.05) (Table II). Ae. sticticus and Ae. vexans had similar 
larval size (6.3 and 6.2 mm; Wilcoxon test, p > 0.05) 
while Ae. rusticus larvae were significantly 70% larger 
than larvae from other species (10.6 mm; Wilcoxon test, 
p < 0.05) (Table II). Correlating larval size with Bti tox- 
ins tolerance indicated that the higher tolerance of Ae. 
rusticus larvae compared to the two other species is due 
to the size differences (Supplementary data), linked to 
its particular life cycle containing a long winter larval 
period. This induces a slow larval development and me- 
tabolism that could lead to an alteration of feeding be- 
haviour, therefore reducing the quantity of Bti ingested 
that could partly explain the lower susceptibility to Bti 
observed for this species. These results are congruent 
with observations made by the local mosquito control 
agency EID that Ae. rusticus requires a higher dose of 
Bti than other Aedes species for efficient control. 

Intraspecific comparisons - Among the populations 
analysed, no increased tolerance to Bti was found in Ae. 
rusticus, Ae. sticticus and Ae. vexans (Figs 2-5). This re- 
sult is congruent with the observations of the EID who 
reported no decreased efficacy of Bti treatments in the 
treated sites over time and also with previous works con- 
ducted in other 5?/-treated areas (Liu et al. 2004, Vasquez 
et al. 2009, Loke et al. 2010, Kamgang et al. 2011). 

The untreated P14 population of ^e. rusticus was sig- 
nificantly 1.77-fold more tolerant than treated P3 and P5 
populations to Cry4Aa, but 1.72-fold less tolerant than 
P5 to CryllAa (Fig. 2). No significant differences be- 



tween these populations were found when tested with 
Cry4Ba and 5?/ (Fig. 2). 

The tolerance of Ae. vexans populations for which 
full bioassays were performed did not show significant 
variations for Cry4Aa, Cry4Ba and Bti while the P44 
population was 2.5-fold more tolerant than the P20 pop- 
ulation to Cryl 1 Aa (Fig. 3). P30 and P42 populations, for 
which the diagnostic doses approach was used, did not 
show any significant differential tolerance to Cry4Ba, 
CryllAa and 5?; (Fig. 4). 

Bioassays with individual Cry toxins being more 
sensitive than bioassays with the whole Bti toxins mix- 
ture, slight changes in tolerance were observed to Cry4- 
Aa and CryllAa in the P14 population of Ae. rusticus 
and to CryllAa in the P44 population of Ae. vexans. 



NS 




Fig. 2: lethal concentrations for 50% individuals (LCj,,) of P03, P05 
and P14 Aedes rusticus populations to Cry4Aa, Cry4Ba, CryllAa 
toxins and commercial Bacillus tfiuringiensis subsp. israelensis (Bti). 
Untreated populations are indicated in white and treated populations 
in blacic. Error: 95% confidence intervals calculated with the probit 
statistical analysis. 1160, 2320 and 2320 larvae were sampled from 
P03, P05 and P14 populations, respectively. Approximate distance be- 
tween the sites: P03-P05, 70 km, P03-P14, 58 km and P05-P14, 85 km. 
Asterisks mean p < 0.05. NS: not significant. 
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P18 1 P20 



Fig. 3: lethal concentrations for 50% individuals (LC,,,) of P18, P20 
and P44 Aedes vexans populations to Cry4Aa, Cry4Ba, CryllAa tox- 
ins and commercial Bacillus tiniringiensis subsp. israelensis (Bti). 
Error: 95% confidence intervals calculated with the probit statisti- 
cal analysis. 380, 2004 and 2000 larvae were sampled from P18, P20 
and P44 populations, respectively. Approximate distance between the 
sites; P18-P20, 21 km, P18-P44, 45 km and P20-P44, 66 km. Asterisks 
mean p < 0.05. NS: not significant. 




Fig. 5: lethal concentrations for 50% individuals (LC^u) of P16, P20 
and P52 Aedes sticticus populations to Cry4Aa, Cry4Ba, CryllAa 
toxins and commercial Bacillus thuringiensis subsp. israelensis (Bti). 
Error: 95% confidence intervals calculated with the probit statistical 
analysis. 1060, 1120 and 2280 larvae were sampled from P16, P20 
and P52 populations, respectively. Approximate distance between the 
sites: P16-P20, 45 km, P16-P52, 58 km and P20-P52, 14 km. Asterisks 
mean p < 0.05. NS: not significant. 




NS 



P30 I P42 
16.3 ng/mL 
C[v4Ba 




NS 



P30 I P42 
32.5 ng/mL 
CryllAa 




P30 I P42 
0 1 jjg/mL 
Bti 



Fig. 4: differential mortality of P30 and P42 Aedes vexans populations 
to Cry4Ba (16.3 ng/mL), CryllAa (32.5 ng/mL) and commercial Ba- 
cillus thuringiensis subsp. israelensis (Bti) (0.1 |ig/mL). Error: stand- 
ard error. 600 larvae were sampled from P30 and P42 populations. 
Approximate distance between the sites: P30-P42, 66 km. NS: not 
significant (Wilcoxon test). 



Even if these changes were significant, tolerance ratios 
were low (less than 2.5-fold) possibly due to natural tol- 
erance variations between populations rather than to an 
increased resistance selected by insecticide treatments. 

The Ae. sticticus population P52 was significantly 
3.5-fold and 8-fold more tolerant to Cry4Aa and Cryll- 
Aa compared to P16 and P20 populations (mean of the 
two populations), while no significant difference was 
found with Bti (Fig. 5). This trend is similar to the re- 
sistance pattern observed for the resistant Aedes aegypti 
strain selected in the laboratory with field-persistent Bti 
sampled in the Rhone -Alpes region which exhibited a 
moderate level of resistance to Bti (3.5-fold), but a higher 
resistance to individual Cry toxins (up to 60-fold) (Tet- 
reau et al. 2012b). It is therefore possible that even if no 
difference of tolerance to Bti was detectable in the three 



populations of ^e. sticticus, a cryptic resistance may be 
developing in the P52 population, only detectable with 
bioassays using each Cry toxins separately. Moreover, 
these fluctuations could be partly explained by the envi- 
ronmental parameters in the collected sites (e.g. pH of the 
water, pollutants etc.) that could affect the susceptibility 
of the larvae, for example by inducing the production 
of detoxification enzymes or metabolites that may inter- 
fere with larval response to Bti toxins (Riaz et al. 2009, 
Poupardin et al. 2012). These effects were limited by 
keeping larvae 24 h in fresh water before performing the 
bioassays. This regularly 5?/-sprayed breeding site was 
located behind a deserted windmill, close to both treated 
and untreated sites. No obvious characteristic, in terms 
of vegetation, ecological corridors and insecticide treat- 
ments strategy could explain the increased tolerance to 
Bti toxins in this particular site. Regular bioassays with 
Bti and individual toxins will confirm if tolerance of this 
population is increasing along with Bti treatments. 

Toward a new diagnostic tool for Bti-resistance de- 
tection - Our results show that the tool developed in the 
present study is more sensitive than standard bioassays 
with Bti. The main limitation of its wide use is the low 
number of larvae sometime available in mosquito breed- 
ing sites. For this reason, the LC^^ determination for each 
Cry toxin, even if it remains the best approach to pre- 
cisely quantify the tolerance ratio between populations, 
does not seem well adapted to highly fragmented small 
and heterogeneous populations. Indeed, the diagnos- 
tic dose approach, which requires fewer toxin amount 
and fewer larvae, seems to be more adapted in order 
to gather tolerance data for a large number of popula- 
tion in a reasonable time. Based on the present study, 
empirical diagnostic doses can be estimated for vari- 
ous species studied and can be used for further studies 
aiming at monitoring Bti resistance in natural mosquito 
populations. Considering that previous works correlated 
a low level of resistance to Bti with high level of resis- 
tance to individual Cry toxins (Paris et al. 2011, Tetreau 
et al. 2012b), it is expected that bioassays with Cry tox- 
ins would detect an incipient resistance to Bti before it 
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could be detected with classical bioassays with Bti. This 
hypothesis has to be confirmed by monitoring the evo- 
lution of the tolerance of the populations tested in this 
study for several years. In combination with bioassays 
with Bti, bioassays with separate Cry toxins allow a 
more sensitive monitoring of the susceptibility of field- 
mosquito populations to Bti. Therefore, we believe that 
such approach should be used in further studies looking 
for 5?/-resistance in other species from the Rhone -Alpes 
region (e.g. Ae. cantans), from the Mediterranean region 
(e.g. Aedes detritus, Aedes caspius) and should also be 
applied to human pathogens vector species such as Ae. 
aegypti, Aedes albopictus and Cx. pipiens exposed to 
regular Bti treatments. 

Bti is widely used for mosquito control. Based on ob- 
servations of the resistance phenotype of a i??/-resistant 
laboratory strain, we developed a new tool: bioassays with 
individual Bti Cry toxins. This highly sensitive tool could 
be used in combination with bioassays with Bti in order 
to improve the monitoring of the susceptibility of field- 
populations of mosquitoes regularly treated with Bti. 
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Lethal concentration for 50% of individuals (LC^^) of Aedes rusticus (circle), Aedes sticticus (square) and Aedes vexans (triangle) populations 
for Bacillus thuringiensis subsp. israelensis (Bti) toxins Cry4Aa (A), Cry4Ba (B), CryllAa (C) and commercial Bti (D) in function of the mean 
larval size (in mm) per population. LC^^ and their 95% confidence intervals (CI) were calculated using a Probit regression. Errors for the larval 
size represent standard deviation (SD) of mean. 



